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ABSTRACT 

Radiance  values  for  the  eight  water  vapor  channels  on  the  DMSP/SSH 
infrared  sounder  were  simulated  for  cloud  free  conditions.  The  simulated 
values  were  inverted  using  a statistical  multiple  regression  procedure 
to  estimate  total  precipitable  water.  Inversions  were  also  obtained 
using  both  the  infrared  water  vapor  channel  data  and  temperature  profile 
data.  For  the  sounder  channels  incorporated  into  the  DMSP/SSH  sensor 
package,  profile  measurements  are  not  possible  below  the  height  of  the 
tropopause  using  the  infrared  channel  data  alone.  Successful  inversions 
with  more  than  an  80  percent  reduction  of  variance  from  climatological 
values  were  accomplished  using  both  temperature  profile  data  (perhaps 
obtained  from  the  CO.,  channel  data)  and  one  or  more  infrared  water  vapor 
channel  observations. 
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1.  INTRODUCTION 


1.1  Objectives 

The  Defense  Meteorological  Satellite  Program  (DMSP)  Block  VD  satel- 
lites will  each  have  a Supplementary  Sensor  H (SSH)  package  with  an 
infrared  sounder  having  channels  in  the  18  to  30  pm  absorption  band  of 
water  vapor.  The  SSH  package  will  also  have  channels  in  the  11  to  15  ym 
absorption  band  of  carbon  dioxide  (CO^)  for  temperature  sounding  and  a 
single  channel  at  10  ym  for  sensing  ozone.  The  objective  of  the  work 
reported  herein,  conducted  by  Environmental  Research  § Technology,  Inc. , 
(ERT) , was  to  develop  an  algorithm  for  the  extraction  of  moisture  para- 
meters from  the  water  vapor  channels  data.  The  moisture  parameters  are 
to  be  used  for  the  correction  of  sea  surface  temperature  maps  and  temper- 
ature profiles  obtained  from  satellite-borne  infrared  radiometers. 

1.2  Summary  of  Results 

The  radiances  for  the  eight  water  vapor  channels  on  the  DMSP/SSH 
package  were  simulated  for  cloud-free  conditions  using  (1)  the  high 
resolution  transmittance  calculations  provided  by  the  Air  Force  Geo- 
physics Laboratory*  (McClatchey  et  al.,  1973);  (2)  the  temperature  and 
water  vapor  profiles  used  for  the  transmittance  calculations;  and 
(3)  the  frequency  dependences  of  each  of  the  channel  filters.  Calcula- 
tions were  performed  for  a set  of  220  atmospheric  profiles  selected  by 
the  Air  Force  Global  Weather  Central  (Wachtmann,  1975)  to  be  statis- 
tically representative  of  the  expected  range  of  temperature  and  moisture 
profiles. 

The  simulated  radiances  were  inverted  using  a statistical  multiple 
regression  procedure  (Gaut  et  al.,  1972)  to  estimate  total  precipitable 
water  and  precipitable  water  profiles  for  the  lower  layers  of  the  atmos- 
phere. Inversions  were  made  using  combinations  of  water  vapor  channel 
and  temperature  profile  data.  Results  for  the  estimation  of  total 
precipitable  water  are  shown  in  Table  1.  The  ratio  of  variances  is  the 

*Data  tapes  with  calculated  transmittance  values  for  220  atmospheric 
profiles  were  provided  by  Dr.  R.  B.  McClatchey  of  AFGL. 
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TABLE  i 

REDUCTION  IN  VARIANCE  FOR  TOTAL  PRECIPITABLE  WATER  ESTIMATORS 


Infrared  Channel 

Number  of 
Temperature  Levels 

Ratio  of  Variances 

8 

S+ 

8.0 

4 ** 

5 

8.1 

7* 

5 

7.2 

8 

0 

2.1 

7* 

0 

1.2 

4** 

0 

2.3 

c--  -1 

5 as  cm 

0 

1.5 

0 

5 

3.2 

C7C  -1 
53s  cm 

5 

7.0 

* Did  not  include  535  cm  1 channel 
**  398,  410,  442  and  535  cm  1 channels 


100,  920,  850,  700,  and  500  mb  temperature  values 


(MViftONMf  NTAl  AfSlAAOft  TCCHHOlOG*  <#C 


a priori  variance  of  the  entire  data  set  divided  by  the  residual  vari- 
ance after  multiple  regression  analysis  or  the  factor  by  which  the 
variance  was  reduced. 

Channel  selection  included  or  excluded  the  1S.~  (535  cm  ) 

channel.  This  channel  is  not  included  on  the  first  SSH  flight  model  but 
will  be  available  on  succeeding  SSH  models.  The  four  channels  at  398, 
110,  442,  and  535  cm  1 were  selected  as  an  optimum  subset  of  the  full 
complement  of  eight  channels.  Finally,  since  the  difference  between 
using  all  eight  channels  and  the  first  seven  channels  was  large,  the  use 
of  only  the  535  cm  * channel  was  explored. 

The  results  listed  in  Table  1 show  that  the  combination  of  eight 
water  vapor  channels  plus  the  temperature  data  at  five  levels  produced 
the  best  results;  the  use  of  the  seven  water  vapor  channels  to  be  flown 
on  the  first  SSH  flight  model  without  including  temperature  data  pro- 
duced no  reduction  in  variance  from  the  a priori  value  obtained  using 
climatological  data.  Results  obtained  using  only  water  vapor  channel 
data  and  including  the  535  cm  1 channel  were  intermediate  between  these 
cases.  The  incorporation  of  the  535  cm  * channel  with  the  temperature 
profile  data  significantly  improved  the  estimate  of  total  precipitable 
water  over  the  use  of  the  single  infrared  channel  alone  or  the  tem- 
perature data  alone. 

The  results  did  not  change  when  both  the  radiance  and  temperature 
values  were  perturbed  by  measurement  noise.  Design  level  noise  of 
0.25  ergs/sec  was  added  to  the  radiance  values  and  a 2°  rms  (root  mean 
square)  temperature  uncertainty  was  used.  The  latter  value  is  typical 
of  the  precision  that  can  be  accomplished  from  infrared  temperature 
profile  measurements  under  cloud-free  conditions. 

The  details  of  the  moisture  parameter  retrieval  analysis  are  pre- 
sented in  this  report.  Section  2 describes  the  inversion  problem  and 
the  statistical  multiple  regression  procedure  used  in  the  analysis; 
Section  3 describes  the  simulation  procedures;  Section  4 provides  a 
discussion  of  results;  and  Section  5 presents  conclusions  and  recom- 
mendations. A description  of  the  software  and  D matrices  for  a sample 
set  of  inversions  are  given  in  the  appendices. 
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2.  THE  INVERSION  PROBLEM 

2.1  The  Radiative  Transfer  Equation 

The  upwelling  thermal  infrared  radiation  sensed  by  radiometers  on 
the  DMSP  satellites  originates  from  surface  emission  and  from  the  atmo- 
sphere above  the  surface.  In  the  absence  of  clouds,  the  radiance 
detected  at  the  satellite  may  be  calculated  using  the  radiative  transfer 
equation  for  an  absorbing  medium  (Chandrasekhar,  1960): 

Iv(o)  = Iv(s)e  Tv(S)  + 0/S  <v(x)  Bv(T(x))e'Tv(x)dx  (2-1) 


where 


Vx) 

is 

radiance  at  x along  the  path  measured  from  the 

satellite 

VX) 

is 

optical  depth  = o/SKy(y)dy 

<v(x) 

is 

absorption  coefficient 

CjV3 

Bv(T(x)) 

= 

r= = Planck  radiance 

(ec2”/T  - 1) 

s 

is 

distance  to  surface 

A 

is 

wavelength 

V 

is 

the  wavenumber 

T 

is 

temperature  (kinetic) 

C1 

= 

2hc  = first  Planck  constant 

C2 

= 

hc/k  = second  Planck  constant 

h 

is 

Planck's  constant 

k 

is 

the  Boltzmann  constant 

Using  the  hydrostatic  equation  and  assuming  that  the  surface  radiates 
as  a black  body,  this  expression  may  be  simplified  to 
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I 

v 


B [T  ]e 

vl  s 


-T  (Ps)  S -T  (P) 

V + J Bv(T(P^e  V 


' K 

V 

PPg, 


dP 


(2-2) 


where 

TV(P) 

P 

P 

P 

S 

T 

s 

g 

V 

<1< 


P k dP 

= / — 
oJ  ppg 

= p(P)  = density 
is  pressure 
is  surface  pressure 

is  surface  temperature 

is  acceleration  due  to  gravity 
is  cos  ii 
is  zenith  angle 


and  the  atmosphere  is  assumed  to  be  horizontally  layered  over  a flat 
surface  (plane  parallel).  Finally,  noting  that 


dr 


dP 

v ppg 


the  radiative  transfer  equation  can  be  reduced  to  the  familiar  form 

Ps  d0  (P) 

I = B [T  ] 0 (P  ) - / B [TCP)  ] —35 — dP  (2-3) 

v vl  sJ  v s o1  vl  J dP 


where 

-T 

0 = e = transmittance 

v 

The  three  forms  of  the  radiative  transfer  equation  were  explicitly 
presented  to  emphasize  the  complexity  of  the  relationship  between  1^ 
measured  at  the  satellite  and  the  atmospheric  parameters  of  interest. 
Moisture  parameters  such  as  mixing  ratio  and  precipitable  water  do  not 
explicitly  appear.  In  the  18  to  30  urn  (330  cm  1 < v < 550  cm  *)  band 
the  only  important  source  of  absorption  is  water  vapor.  The  absorption 
coefficient  depends  upon  the  contributions  of  a number  of  rotational- 
vibrational  absorption  lines  in  the  vicinity  of  v.  Following  Burch 


6 


r 


— — — - 


mvmowiia'x  KisioCHt  >irB»oioc* 


et  al  (1974),  the  absorption  coefficient  may  be  modeled  as  a sum  of  the 
contribution  of  nearby  lines  (local)  and  a continuum  contribution  contri- 
buted by  the  far  wings  of  a large  number  of  distant  lines.  The  local 
pressure  broadened  lines  may  be  described  by  the  classical  Lorentc  line 
shape  and  the  contribution  of  each  line  may  be  calculated  given  the 
frequency,  width  and  intensity  of  each  line.  The  Lorent;  line  shape 
does  not  adequately  represent  the  actual  line  intensity  more  than  10  - 15  cm 
from  the  line  center.  At  larger  frequency  offsets,  the  contributions 
from  a large  number  of  distant  lines  becomes  important  and  their  contri- 
butions may  be  modeled  as  background  continuum  absorption. 

The  model  for  the  absorption  coefficient  is  given  by: 


KT 

v 


c 


e 

pkT 


[Z  C.p  + C e + C p] 
. i*  s n*  1 


(2-4) 


where 

e is  the  partial  pressure  of  water  vapor, 

is  the  continuum  coefficient  for  self  broadening, 

is  the  continuum  coefficient  for  total  pressure  (nitrogen) 
broadening,  and 

C.  are  the  contributions  by  lines  at  v.  with  the  line  width 

l 7 l 

pressure  dependence  made  explicit. 

Noting  that  e/T  <*  p , the  water  vapor  density,  Equation  2-4  may  be 
represented  by 

K 

^ = r [f 1 (v,P,T)  ♦ r f 2 (v , P ,T) ] 

to  make  the  dependence  on  mixing  ratio,  r,  explicit.  Now  Equation  2-2 
may  be  rewritten 

-T  (P  ) PS  -T  (P) 

I =B(T)eVS  + / B [T ( P ) ] e V r(f.  ♦ rf,)—  (2-5) 

v v s oJ  v 1 ‘ 1 2 ug 
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where 


,P  tc  c ,dP 
J r(f.  + rf.) — 
o 1 2 ug 


The  dependence  of  the  measured  radiance  on  mixing  ratio  is  highly 
nonlinear.  The  optical  depth  is  a temperature  and  pressure  weighted 
integral  of  the  mixing  ratio  plus  a term  with  a second  order  dependence 
on  mixing  ratio.  The  radiance  values  have  a further  nonlinear  dependence 
on  the  weighted  integral  of  mixing  ratio  through  the  exponential  func- 
tion. Equation  2-5  is  a nonlinear  equation  for  r when  I is  given.  To 
obtain  information  about  r or  about  precipitable  water,  W, 


W 


(P) 


p 

r rdP 
oJ  g 


the  equation  must  be  inverted. 


2.2  Statistical  Parameter  Inversion  Method 


The  integral  equation  (Equation  2-5)  is  highly  nonlinear  in  r. 

Smith  and  Howell  (1971)  attempted  to  invert  this  equation  using  an 
iterative  technique.  To  be  sucessful,  the  iterative  solution  should 
take  into  account  the  variations  of  the  weighting  functions  with  the 
mixing  ratio  profile  and  temperature  profile.  For  the  analysis  of 
temperature  profiles  using  iterative  techniques  such  as  the  minimum 
information  technique,  the  variation  of  the  weighting  functions  with  the 
parameter  of  interest,  temperature,  is  assumed  to  be  small  and  a fixed 
weighting  function  can  be  assumed.  This  is  not  possible  for  the  inver- 
sion of  water  vapor  data.  The  result  is  a lengthy,  time  consuming  pro- 
cess which,  in  the  end,  may  not  converge  to  the  desired  result  due  to 
the  non-uniqueness  of  the  solution. 

The  Statistical  Parameter  Inversion  Method  developed  by  ERT 
(Gaut  et  al.,  1972)  bypasses  the  requirement  for  the  lengthy  iterative 
procedure  for  the  inversion  of  the  radiative  transfer  equation  by  using 
multiple  regression  techniques.  Multiple  regression  procedure  has  been 
used  by  many  investigators  for  the  inversion  of  the  radiative  transfer 
equation.  Recently,  Smith  and  Woolf  (1976)  described  a procedure  nearly 
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identical  to  the  one  used  by  Gaut  et  al.  (1972)  for  the  inversion  of 
Nimbus-6  data  to  generate  temperature  and  mixing  ratio  profiles.  As 
applied,  a multiple  regression  analysis  is  performed  using  radiance  data 
or,  in  their  absence,  simulated  radiances  as  the  independent  set  of 
variates  and  the  desired  parameters  as  the  dependent  set.  A linear 
relationship  is  assumed  between  the  dependent  and  independent  sets  and 
the  least  square  and  best  estimate  coefficients  relating  the  two  sets 
are  determined  which  minimize  the  differences  between  the  estimated  and 
observed  values.  This  procedure  selects  the  most  probable  set  of  para- 
meters from  the  infinite  set  of  parameters  that  could  have  produced  the 
same  radiance  values. 

Let  v^  be  the  ith  parameter,  say  W^,  and  d^  be  the  j1*1  element  of 
a measured  data  set,  then  the  assumed  relationship  is 


v : = 


N 
E D. 


d. 


j-1  1J  3 


which,  in  matrix  notation,  is 


v = £ d 


where 

v7  is  the  estimated  parameter, 

D.j  are  the  elements  of  a linear  operator  n relating  d^  to  v7,  and 
N is  the  number  of  data  elements  in  an  observation. 


The  problem  is  to  determine  the  elements  of  the  linear  operator  such 

that  the  variance  of  the  difference  between  v.  and  v*  is  a minimum 

1 i 

(least  square  minimum  error).  The  elements  of  the  linear  operator  can 
be  estimated  using  multiple  regression  techniques  and  simultaneous  sets 
of  observations  (or  simulations)  of  v and  d.  The  elements  £ then  minimize 
the  error  for  the  sets  of  observations.  For  use  as  an  estimation  proce- 
dure with  a new  set  of  observations,  the  elements  of  £ will  only  provide 
minimum  variance  if  the  statistical  properties  of  both  the  original  and 
new  sets  are  identical.  This  is  usually  taken  on  faith  if  simultaneous 
v and  d elements  are  not  available  or  can  be  tested  if  they  are  available. 
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The  Statistical  Parameter  Estimation  Method  first  determines  the 
empirical  orthogonal  functions  that  best  represent  the  data  elements. 

The  data  are  then  represented  by  a truncated  set  of  the  empirical  ortho- 
gonal functions  to  minimize  the  effect  of  noise  and  redundant  observa- 
tions. The  latter  set  is  used  in  a classical  regression  analysis  to 
determine  the  D matrix  elements.  The  calculations  are  made  using  a 
modified  d'  vector  with  the  first  element  set  to  unity  and  the  remaining 
elements  equal  to  the  departures  from  their  expected  mean  values  nor- 
malized by  the  expected  mean  value 


{d!  > = 

x 


^i-l" 


<di-l> 


where  the  < > operation  represents  averaging  over  the  number  of  observa- 
tions in  the  set  used  to  generate  the  D matrix.  The  correlation  matrix, 

Q (d',d')  where  . = <d|  dj>,  generated  using  the  values  for  d'  and 

is  diagonalized  to  determine  both  its  eigenvalues  and  eigenvectors.  The 
eigenvectors  are  the  empirical  orthogonal  functions  for  the  data  elements 
and  the  eigenvalues  give  the  amount  of  the  variance  represented  or 
explained  by  each  eigenvector  (Lorentz,  1956). 

T 

The  transpose  of  the  eigenvector  matrix  (R  ) is  used  to  transform 

T 

the  modified  data  vector  d'  to  a new  data  vector  d"  = R d'  which  may  be 
truncated  to  provide  a best  representation  of  the  data  vector  for  use  in 
the  regression  analysis.  Using  the  truncated  set. 


D"  = C(v, d")  £'l(d",d") 


The  H'  matrix  for  use  in  estimating  v is  then  R'  = R"*R  . The  R 
matrix  calculated  in  this  manner  minimizes  sensitivity  to  noise  and 
correlation  between  data  elements;  the  correlated  data  elements  being 
represented  by  a single  component  of  the  empirical  orthogonal  function 
set  ((laut  et  al,  19~2;  Smith  and  Woolf,  1976). 
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3.  SIMULATION  OF  RADIANCES 

5.1  Band  Pass  Filters 

The  use  of  the  ()'  matrix  for  the  inversion  of  infrared  sounder 
measurements  first  requires  the  availability  of  both  data  and  simul- 
taneously observed  parameters  for  the  generation  of  a matrix.  In  the 
absence  of  actual  observations,  simulated  data  must  be  used.  The  radi- 
ances to  be  observed  in  each  of  the  channels  of  the  DMSP/SSH  sounder  are 
not  single  frequency  observations  as  implied  in  the  equations  of  Section  2.1 
but  are  spaced  over  the  frequency  band  defined  by  the  instrument  filters. 

Table  2 describes  the  frequency  bandpass  characteristics  of  the 
F channels  of  the  first  flight  model  water  vapor  sounding  (Barnes 
Engineering  Co.,  1975)  and  of  the  eighth  channel  to  be  included  on  the 
subsequent  flight  models.  The  numbers  1 to  8,  given  in  Table  2,  are 
used  to  refer  to  the  channels.  Their  bandpass  characteristics  in  terms 
of  energy  throughput  or  transmission  are  displayed  on  Figure  1.  The 
data  on  this  figure  were  compiled  from  measurements  at  0.5  cm  1 spacings 
provided  by  the  manufacturer  (Barnes  Engineering  Co.,  1975)  and  were 
used  to  define  the  response  of  each  of  the  sounder  channels.  Although 
the  full  system  frequency  response  as  well  as  the  filter  response  were 
measured  by  the  manufacturer,  digitized  data  were  available  for  the 
filters  and  were  used  for  the  simulation.  As  shown  in  Table  2,  the 
differences  between  using  the  two  sets  of  data,  as  calculated  by  the 
manufacturer  and  as  calculated  from  the  digitized  filter  data  plotted  on 
Figure  1 are  small.  The  data  for  channel  8 were  obtained  from  measure- 
ments on  a similar  channel  filter  on  the  sensor  on  satellite  DMSP-6530. 

The  radiances  to  be  measured  on  the  satellite  are  related  to  the 
theoretical  values  described  in  Section  2.1  by 

oo 

I . = / g . (v) I dv 

l ‘ &i  v 

o 

where  g.  is  the  transmission  function  for  the  i*^  channel. 
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534.35  11.55  18.43 


• Ill  42 


Figure  1 Wavenumber  Dependence  of  the  SSH  Instrument  Filters 
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3.2  Calculations  of  Transmittance 


The  transmittance  values  needed  at  each  frequency  for  the  calcula- 
tion of  I.  require  the  calculation  of  the  absorption  coefficient  for 
each  of  the  lines  within  the  interval  and  for  the  continuum  contribu- 
tions from  distant  lines.  The  transmittance  values  for  each  line 
0^(p)  = e VP)  were  provided  by  Dr.  R.  McClatchey  of  AFGL  for  0.1  cm  * 
increments  and  12  pressure  levels  (McClatchey  et  al.,  1973).  The  cal- 
culations were  made  for  220  atmospheric  profiles  selected  by  R.  F. 
Wachtmann  of  Air  Force  Global  Weather  Center  (AFGWC)  to  be  statistically 
representative  of  the  profiles  to  be  expected  on  an  operational  basis. 

He  compiled  a larger  set  of  profiles  and  determined  that  this  subset  was 
statistically  identical  with  the  larger  set  (Wachtmann,  1975).  The  set, 
therefore,  is  useful  for  the  estimation  of  a £ matrix  because  the 
statistical  representativeness  of  the  set  was  considered  in  its  selection. 

The  calculations  provided  by  Dr.  McClatchey  did  not  include  the 
continuum  contributions.  These  were  established  from  observations  made 
by  Burch  et  al.  (1974).  The  observations  were  made  at  two  temperatures 
in  spectral  gaps  between  337  and  822  cm  * . The  observations  for  Cg  and 
Cn  are  presented  in  Figures  2 and  3.  The  data  were  fit  by  smooth  curves 
that  matched  the  lowest  values  at  the  reported  frequencies.  This  pro- 
cedure was  followed  because  the  data  could  be  contaminated  by  the  edges 
of  neighboring  lines.  The  minimum  values  should  represent  the  continuum 
contribution  with  the  least  amount  of  contamination  by  local  effects. 

The  curves  were  separately  established  for  each  temperature,  and  values 
were  scaled  from  each  curve  for  the  center  frequency  of  each  channel  and 
tabulated  (Table  3).  Noting  that  line  broadening  can  be  modeled  as 

- 

proportional  to  Tn  (Goody,  1964),  the  two  temperature  values  were  used 
to  establish  a value  of  n that  best  fits  all  the  data  for  each  coeffi- 
cient. For  self-broadening,  n = -2;  for  nitrogen  broadening,  n = *2. 

These  data  were  used  to  calculate  the  continuum  contribution  to  the 

i 

transmittance.  From  Equation  2-3,  it  is  evident  that  the  total  trans- 
mittance is  given  by 
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TABLE  3 

WATER  VAPOR  CONTINUOUS  ABSORPTION  COEFFICIENTS 


Channel 

Nominal 

Center 

(cm*1) 

Optical 

Center 

(cm-1) 

Absorption  Coefficient 

C*. 

(430  K) 

Cs 

(296#K) 

Cn 

(430*K) 

cn 

(296°K) 

1 

3S3.0 

353.78 

5757 

12150 

1245 

590 

2 

359.0 

356.67 

5450 

11400 

1055 

500 

3 

373.5 

373.46 

4194 

8850 

654 

310 

4 

397.  S 

398.15 

3270 

6900 

369 

175 

5 

409.0 

409.83 

2962 

6250 

285 

135 

6 

423.0 

418.96 

2725 

5750 

232 

110 

7 

441.0 

442.19 

2227 

4700 

137 

65 

8* 

534.5 

534.35 

972 

2050 

53 

25 

•Not  on  Flight  I/SSH  (Filter  Response  Function  from  DMSP-6530) 


where  0.  is  the  transmittance  for  the  lines  within  the  filter  bandpass 

1 , X, 

weighted  by  the  filter  response 


0i,£(P) 


' 

= /Qv(P)gi(v)dv 


Jg, (v)dv 


and  0.  (P)  is  the  continuum  transmittance  evaluated  at  the  center 

i » c 

frequency  for  the  channel.  The  filter  weighted  line  transmittance 
values  0^  ^(P)  were  calculated  using  the  0^  ^(P)  values  at  0.1  cm 
intervals  and  by  a second  order  spline  interpolation  between  the  0.S  cm 
filter  response  values. 

The  continuum  calculations  were  made  for  each  channel  center  fre- 

2/7 

quency  using  35  pressure  levels  equally  spaced  in  units  of  P between 
0.01  and  1,000  mb.  The  values  of  T,  P and  r required  to  make  the  cal- 
culations were  the  average  values  between  the  layers.  The  transmittance 
values  for  the  lines  were  also  spline  interpolated  to  the  35  pressure 
levels.  The  total  transmittance  values  were  calculated  for  each  channel 
and  pressure  level  for  use  in  calculating  the  simulated  radiances. 
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Weighting  functions  were  calculated  for  several  of  the  profiles  to 
display  the  results  of  the  transmittance  calculations.  The  weighting 
functions  for  the  Planck  radiance  (Equations  2-1  and  2-3)  were  cal- 
culated by  the  weighting  function 

d0.  d0. 

1 l 

dH  d(lnP) 

Weighting  functions  for  all  the  channels  for  a tropical  summer  atmos- 
phere are  displayed  on  Figure  4.  With  the  exception  of  channel  8,  all 
the  weighting  functions  peak  at  500  mb  or  above.  Although  these  weight- 
ing functions  are  not  the  coefficients  of  precipitable  water  or  of 
mixing  ratio,  the  parameters  to  be  sensed  (Smith  and  Howell,  1971),  they 
show  that  only  the  Planck  radiances  in  the  upper  layers  of  the  atmos- 
phere contribute  to  the  radiance  values.  As  displayed,  the  weighting 
functions  contain  the  water  vapor  information.  These  functions  show 
that  only  water  vapor  above  500  mb  can  be  sensed  for  this  profile. 

A subset  of  four  channels  were  selected  that  had  weighting  functions 

that  peaked  at  different  heights  in  the  lower  atmosphere.  Weighting 

functions  for  channels  4,  5,  7 and  8 at  398,  410,  442  and  534  cm  1 , 

respectively,  are  displayed  for  two  widely  different  atmospheres  on 

Figures  5 and  6.  Figure  5 is  for  the  same  atmospheric  profile  as  Figure  4, 

tropical  summer.  Figure  6 represents  an  arctic  winter  profile.  These 

data  show  that  channels  4,  5 and  7 peak  at  the  same  height  for  the 

moist,  warm  tropical  summer  profile.  The  weighting  functions  peak  at 

much  lower  heights  for  the  cold,  dry  arctic  winter  profile.  In  the 

latter  case,  channel  8 should  provide  useful  information  down  to  the 

surface.  These  data  display  the  wide  variability  possible  for  the 
_ t m 

product  x^e  v (see  Equation  2-1). 
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Pressure  (mb) 


Figure  4 


Pressure  Dependence  of  the  Weighting  Functions  for  the 
Planck  Radiance  for  Each  Water  Vapor  Channel 
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3.3  Radiance  Calculations 

The  radiance  values,  I , were  calculated  using  the  35  pressure 
levels  and  evaluating  the  Planck  raGiances  for  the  layers  between  the 
radiance  levels.  From  Equation  2-3,  it  is  evident  that  a surface 
temperature  must  be  specified  to  calculate  the  surface  emission  for  use 
in  the  radiative  transfer  equation.  Surface  temperature  measurements 
were  not  provided  with  the  220  atmospheric  profiles.  Surface  tempera- 
ture values  for  use  in  the  equation  were  randomly  generated  to  be  within 
1.5°  rms  of  the  1,000  mb  temperature  value.  The  surface  is  assumed  to 
be  at  1,000  mb  and  to  have  unity  emissivity. 

The  radiance  calculations  were  made  for  each  of  the  220  profiles. 
Calculated  radiance  values  for  the  two  samples,  the  arctic  winter  and 
tropical  summer  profiles,  are  displayed  on  Figure  7.  These  data  show 
large  variations  between  the  observations.  The  data  are  of  the  same 
order  of  magnitude  as  the  arctic  and  tropical  observations  reported  by 
Hanel  and  Conrath  (1970). 
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Figure  7 Radiance  Estimate  Summary  for  the  220  Atmospheric  Profiles 
and  Examples  for  Arctic  Winter  and  Tropic  Summer 
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4.  ANALYSIS 


4.1  Sample  Variability 

The  set  of  220  simulated  radiances  and  their  associated  atmospheric 
profiles  comprise  the  data  and  parameter  set  used  for  the  generation  of 
the  |)  matrix.  The  utility  of  different  retrieval  algorithms  using 
different  channels  and  combinations  of  channel  radiance  values  and  tem- 
perature values  is  measured  by  the  factor  by  which  the  a priori  or 
climatological  variance  of  the  parameter  is  reduced.  The  variance 
values  and  covariances  required  to  generate  the  £ matrices  used  for  the 
retrieval  of  precipitable  water  were  calculated  using  the  entire  set  of 
220  atmospheric  profiles.  The  profiles  were  selected  from  a larger  set 
and  constrained  to  have  the  same  statistical  properties  as  the  larger 
set  (Wachtmann,  1975).  The  resultant  data  set  represents  a number  of 
profiles  from  different  latitude  regions  and  months  of  the  year  as  shown 
in  Table  4.  The  mean  and  standard  deviation  of  the  calculated  radiances 
for  this  data  set  are  shown  on  Figure  7. 


TABLE  4 

DISTRIBUTION  OF  SOUNDINGS  BY  LATITUDE  AND  NORTHERN  HEMISPHERE  MONTHS 


Latitude 

Month 

1-1S0 

15-30* 

30-45° 

45-60* 

60-75* 

75-90* 

Jan 

16(10) 

18 

14 

37 

20 

13 

Feb 

10(4) 

18 

17 

25 

10 

15 

Mar 

14(4) 

20 

17 

18 

9 

10 

Apr 

8(8) 

17 

12(2) 

29 

18 

14 

May 

18(10) 

21 

24(2) 

12 

15 

16 

Jun 

16(8) 

26 

21(4) 

24 

22 

10 

Jul 

17(11) 

23 

24(2) 

23(3) 

15(9) 

15(4) 

Aug 

15(7) 

20 

20 

23(7) 

25(5) 

12 

Sep 

10(10) 

12 

15 

0 

6(4) 

4 

Oct 

20(9) 

16 

20 

23 

10 

10 

Nov 

13(7) 

18 

16 

21 

10 

10 

Dec 

19(10) 

16 

16 

44 

10 

6 

•The  numbers  in  parentheses  represent  that  part  of  the  total  taken  in 
the  southern  hemisphere.  Southern  hemisphere  soundings  are  counted  for 
the  month  six  months  out  of  phase  of  the  observation  month  in  order  to 
represent  seasonal  distribution. 


# PRECgDIhC  rA££  BLANK- NOT  FILMED 
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The  moisture  parameter  selected  to  test  the  inversion  algorithms 
for  the  simulated  data  are  the  precipitable  water  values  for  five 
pressure  levels  between  1,000  and  500  mb.  The  precipitable  water  is 
defined  by  Equation  2-6.  The  total  precipitable  water  is  the  value  for 
1,000  mb  pressure.  This  value  is  of  most  interest  for  sea  surface  tem- 
perature mapping  and  was  used  by  Cogan  and  Willand  (1975)  as  an  indepen- 
dent variable  in  the  procedure  for  estimating  the  temperature  deficit 
for  the  correction  of  the  sea  surface  values.  Similarly,  the  precipi- 
table water  values  for  different  pressure  levels  are  required  to  esti- 
mate corrections  for  the  transmission  functions  used  in  the  retrieval  of 
temperature  profile  data  from  the  C0o  sounder  channels'  data.  The  mean 
and  standard  deviations  of  precipitable  water  are  depicted  on  Figure  8 
and  are  listed  in  Table  5.  From  the  table,  it  is  evident  that  the 
a priori  or  climatological  standard  deviation  of  precipitable  water 
varies  from  60  to  80  percent  of  the  mean  depending  upon  the  pressure 
level.  The  individual  tropical  summer  and  arctic  winter  profiles  used 
in  the  radiance  calculations  depicted  in  Figure  7 are  displayed  in 
Figure  8. 

The  precipitable  water  values  are  not  all  independent  parameters. 

The  first  three  eigenvectors  in  an  empirical  orthogonal  function  repre- 
sentation of  the  parameter  set  explain  99.91  percent  of  the  variance. 

A similar  analysis  for  the  eight  water  vapor  channels  shows  that  the 
radiances  are  highly  correlated.  The  first  three  eigenvectors  in  an 
empirical  orthogonal  function  representation  of  the  radiances  explain 
99.966  percent  and  the  first  four  eigenvectors  explain  99.997  percent. 

The  temperature  values  at  the  same  pressure  levels  used  for  precipitable 
water  have  the  mean  and  variance  values  listed  in  Table  6.  Again,  an 
empirical  orthogonal  function  representation  of  the  temperature  data 
shows  that  only  two  eigenvectors  are  required  to  explain  99.9'  percent 
of  the  variance.  These  results  show  that  both  the  data  elements  arc 
highly  correlated  and  the  parameter  elements  arc  also  highly  correlated. 
This  means  that  the  number  of  independent  pieces  of  information  in  an 
observation  is  significantly  smaller  than  the  number  of  sounder  channels. 
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Pressure 

Mean 

Variance 

Standard  Deviation 

Standard  Deviation 

(mb) 

2 

(gm/cm~) 

(gm/cm  ) 

(gm/cm‘) 

Percent  of  Mean 

1000 

2.9 

2.7 

1.6 

57 

920 

2.0 

1.3 

1.2 

59 

8S0 

1.5 

0.78 

0.88 

61 

700 

0.66 

0. 19 

0.44 

66 

500 

0. 14 

0.012 

0.11 

77 

Maan  i Standard  Daviation 
O Arctic  W intar 
A Tropic  Summar 


Prccipitdblc  Wcttr  (9/cm2) 


Figure  8 Precipitable  Water  Estimate  Summary  for  the 
220  Atmospheric  Profiles  and  Examples  for 
Arctic  Winter  and  Tropic  Summer 
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TABLE  6 

TEMPERATURE  VALUES 


Pressure 

(mb) 

Mean 

(°K) 

Variance 

(°K)‘ 

Standard  Deviation 
(°K) 

1000 

287.7 

175 

13.2 

920 

285.0 

145 

12.1 

850 

282.3 

128 

11.3 

700 

275.1 

92 

9.6 

500 

259.5 

85 

_ 

9.2 

4.2  Results  of  Precipitable  Water  Retrieval 

The  Statistical  Parameter  Estimation  Method  was  applied  to  a number 
of  combinations  of  sounder  channels  and  temperature  levels  to  determine 
the  optimum  combination  for  use  with  the  DMSP/SSH  sounder.  The  empirical 
orthogonal  function  representation  of  the  data  elements  were  truncated 
to  eliminate  all  eigenvectors  that  explained  less  than  0.01  percent  of 
the  variance.  This  procedure  retained  the  smallest  number  of  elements 
required  to  adequately  describe  the  data  while  reducing  the  effects  of 
measurement  noise.  The  truncation  criteria  was  selected  to  eliminate 
the  effects  of  noise  while  retaining  only  the  significant  data.  A com- 
parison of  results  using  several  truncation  criteria  are  displayed  in 
Table  7.  Using  a truncation  level  of  10  ^ (0.1  percent),  the  variance 

for  total  precipitable  water  was  reduced  only  30  percent.  Reducing  the 

-4 

truncation  level  to  10  (0.01  percent)  improved  the  estimate  of  total 

precipitable  water.  In  this  case,  the  variance  was  reduced  by  53  per- 
cent both  with  and  without  contamination  by  measurement  noise.  The 
measurement  noise  was  used  to  randomly  perturb  the  data  elements  prior 
to  the  generation  of  the  Q matrix.  The  percent  variance  represented  by 
the  noise  of  0.002,  which  is  rejected  by  the  selected  truncation  level. 

A lower  truncation  level  would  produce  noise  sensitive  results. 
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A number  of  sounder  channel  and  temperature  level  parameters  were 
used  in  the  analyses  of  water  vapor  retrievals.  Temperature  was 
included  as  a data  element  because  it  determines  the  maximum  value  of 
precipitable  water  that  the  atmosphere  can  carry.  It  is  often  noted 
that  relative  humidity  values  vary  over  a narrow  range.  The  mixing 
ratio  values,  therefore,  should  be  highly  correlated  with  the  tempera- 
ture values.  Since  the  precipitable  water  is  dominated  by  the  contribu- 
tions of  the  mixing  ratio  at  the  lowest  level  used  in  the  calculation, 

| the  layer  temperature  value  and  precipitable  water  value  should  be 

highly  correlated.  This  was  true  for  the  220  sample  profile  set  as 
indicated  in  Table  8.  The  correlations  between  temperature  and  precipi- 
table water  ranged  from  0.5  to  0.8.  The  correlations  between  precipi- 
table water  and  radiance  values  were  smaller,  ranging  from  0.0  to  0.6. 

A subset  of  four  infrared  sounder  channels  were  also  used  for  the 
estimation  of  precipitable  water.  Channels  4,  5,  7 and  8 (398,  410,  441 
and  535  cm"*)  were  selected  as  having  the  highest  correlation  with 
precipitable  water  (Table  8)  and  the  largest  differences  in  weighting 
functions.  A subset  of  the  original  data  set  was  selected  as  a method 
for  operationally  reducing  the  complexity  of  precipitable  water  retrieval 
algorithms. 

The  results  of  Statistical  Parameter  Estimation  Method  retrievals 
for  each  of  the  channel -temperature  level  combinations  are  presented  in 
Tables  9 through  13.  The  data  on  each  table  are  for  a single  precipi- 
table water  parameter.  The  results  show  that  reductions  in  variance  of 
80  to  90  percent  are  possible  using  temperature  data  at  four  or  five 
levels  and  at  least  one  infrared  sounder  channel.  The  differences 

between  the  use  of  one  or  more  sounder  channels  are  small.  The  use  of 

l 4 

temperature  only  data  provides  a 70  percent  reduction  in  variance  for 
total  precipitable  water  but  using  temperature  in  the  lower  levels  only 
is  not  as  useful  for  precipitable  water  for  levels  above  700  mb. 

Sounder  data  only  is  most  useful  when  channel  8 is  included  but  no 
reduction  in  variance  (at  the  0.05  significance  level)  occurs  when 
channel  8 is  not  included. 
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TABLE  9 


TOTAL  PRECIPITABLE  WATER  RETRIEVAL  W 


(1000) 


Channel 

— 

Temperature 

Levels 

(mb) 

. 

Standard  Deviation 
of  Residuals 
(gm/cra2) 

Ratio  of 
Variance 

Percent  Reduction 
in  Variance 

(%) 

1-8 

1000-500 

0.58 

8.0 

87 

1-8 

1000-700 

0.61 

7.1 

86 

1-7 

1000-500 

0.61 

7.2 

86 

4,5,7  6 8 

1000-500 

0.S7 

8.1 

88 

8 

1000-500 

0.62 

7.0 

86 

No 

1000-500 

0.92 

3.2 

69 

1-8 

No 

1.1 

2.1 

53 

1-7 

No 

1.1 

4,5,7  6 8 

No 

1.1 

2.3 

S6 

4,5  6 7 

No 

l.S 

1.2 

17 

8 

No 

1.3 

L_±i 

35 

TABLE  10 

920  mb  PRECIPITABLE  WATER  RETRIEVAL  W(92Q) 


’ Channel 

Temperature 

Levels 

(mb) 

Standard  Deviation 
of  Residuals 
(gm/cm2) 

Ratio  of 
Variance 

Percent  Reduction 
in  Variance 
(') 

1-8 

1000-500 

0.40 

8.5 

88 

1-8 

1000-700 

0.42 

7.5 

87 

1-7 

1000-500 

0.43 

7.2 

86 

4 , S , 7 6 8 

1000-500 

0.39 

8.6 

88 

8 

1000-500 

0.43 

7.3 

86 

No 

1000-500 

0.73 

2.6 

61 

1-8 

No 

0.84 

1.9 

47 

1-7 

No 

1.1 

1.1 

11 

4,5,7  6 8 

No 

0.82 

2.0 

50 

4,5  6 7 

No 

1.1 

1.1 

11 

8 

No 
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Temperature 

Levels 

(mb) 


1-8 

1-8 

1-7 

4,5,7  6 8 
8 

No 

1-8 

1-7 

4,5,7  6 8 
4,5  6 7 
8 


.28 

9.8 

90 

.30 

8.5 

88 

.32 

7.7 

87 

.28 

10.0 

90 

.31 

8.2 

88 

.59 

2.3 

56 

TABLE  12 

700  mb  PRECIPITABLE  WATER  RETRIEVAL  W 


Temperature 

Standard  Deviation 

Ratio  of 

Percent  Reduction 

Levels 

of  Residuals 

Variance 

in  Variance 

(mb) 

(gm/cm2) 

(%) 

13 

.3 

10. 

.6 

8. 

.8 

13. 

,6 

11. 

.0 

1. 

8 

1. 

.6 

1. 

.0 

1. 

7 

1. 

.0 

1 

,i 

1 8 

1000-500 

1-8 

10O0-700 

1-7 

1000-500 

4,5.7  6 8 1 

1000-500 

8 

1000-500 

No 

1000-500 

1-8 

No 

1-7 

No 

4,5,7  8 8 

No 
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4.3  Applications 


The  water  vapor  retrieval  algorithm  using  both  infrared  and  tempera- 
ture profile  data  is  capable  of  estimating  total  precipitable  water  with 

-> 

an  rms  error  of  less  than  0.6  gm/cnT.  Shen  and  Smith  (1972)  estimated 

that  the  uncertainty  in  the  estimation  of  total  precipitable  water 

2 

achieval  with  the  SIRS-B  sounder  was  0.6  gm/cm-,  the  same  as  estimated 
for  the  DMSP/SSH  sounder  using  the  established  regression  procedure  and 
both  sounder  and  temperature  data.  The  inversions  may  be  made  using  the 
matrix  operation  V = £ d where  the  matrix  elements  are  listed  in 
Appendix  2.  Since  the  220  element  data  set  was  generated  using  statisti- 
cally representative  samples,  the  matrices  generated  using  the  simu- 
lated radiance  values  are  applicable  for  operational  use.  The  £) 
matrices  generated  in  this  fashion  may  not  be  the  best  for  operational 
use  in  that  the  radiances  are  simulated,  not  measured,  and  possible 
errors  in  the  model  can  affect  the  results.  In  practice,  |)  matrices 
generated  using  simultaneous  cloud-free  satellite  observations  and 
profile  measurements  are  to  be  preferred. 

The  water  vapor  correction  technique,  described  by  Cogan  and 
Willand  (1975),  relies  upon  the  use  of  the  total  precipitable  water  as  a 
predictor  of  temperature  deficit.  For  a ^ ^ j 000 ) estimation  with  a 
0.6  gm/cm^  rms  error,  the  uncertainty  in  the  prediction  of  the  tempera- 
ture deficit  is  less  than  0.6°K.  From  the  analysis  presented  by  Cogan 
and  Willand,  the  uncertainty  in  the  estimation  of  sea  surface  temperature 
would  be  less  than  1°K  rms. 

The  procedure  used  for  the  extraction  of  precipitable  water  is 

applicable  to  cloud-free  radiance  observations.  Temperature  profile 

estimates  with  rms  errors  of  less  than  2°K  are  required  for  making 

2 

estimates  with  less  than  0.6  gm/cm  rms  error  in  These  may  be 

obtained  from  the  C0-,  temperature  sounder  for  cloud-free  conditions. 

The  required  data,  therefore,  may  be  entirely  obtained  from  the  DMSP/SSH 
sensors.  The  procedure  applies  only  for  cloud-free  conditions,  hence 
the  presence  of  clouds  must  be  sensed.  Since  sea  surface  temperature 
measurements  can  only  be  made  under  cloud-free  conditions,  this  restric- 
tion on  the  retrieval  of  water  vapor  information  is  not  critical. 
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5.  CONCLUSIONS 

The  220  simulated  radiances  sets  were  used  to  estimate  the  coefficients 
of  a g matrix  linear  operator  for  use  in  estimating  precipitable  water 
using  infrared  radiance  and  temperature  level  data.  The  residual  errors 
expected  using  the  statistical  regression  procedure  are  less  than 
20  percent  of  the  mean  values  for  precipitable  water  in  the  lower  levels 
of  the  atmosphere.  Smith  and  Woolf  (1976)  report  better  results  for 
Nimbus-6  data  when  the  infrared  and  microwave  data  are  combined.  The 
microwave  channels  are  relatively  transparent  and  water  vapor  may  be 
sensed  to  the  surface.  The  weighting  functions  displayed  in  Figures  5 
and  6 show  that  under  warm,  humid  conditions,  water  vapor  near  the 
surface  cannot  be  sensed  by  the  DMSP/SSH  sounder.  The  success  of  the 
combined  sounding  plus  temperature  profile  algorithm  for  warm,  humid 
profiles  is  a result  of  the  natural  correlation  between  temperature  and 
water  vapor.  Improved  water  vapor  retrievals  could  be  obtained  if 
additional  infrared  sounding  channels  were  provided  in  a less  opaque 
region  of  the  water  vapor  spectrum.  A number  of  the  current  F-channels 
on  the  SSH  package  are  redundant  and  may  be  eliminated  without  loss  in 
sensor  capability. 

The  Statistical  Parameter  Estimation  Procedure  was  used  for  water 
vapor  retrievals.  It  can  also  be  successfully  used  for  the  extraction 
of  temperature  profile  data  from  the  C0-,  sounder  channels.  The  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  operational  procedures  are 
now  being  changed  to  use  the  statistical  parameter  estimation  procedure 
for  profile  estimation  (Smith  and  Koolf,  19“6).  \ similar  move  for  the 

extraction  of  data  from  the  DMSP/SSH  data  is  recommended.  The  simultan- 
eous inversion  of  K,  F and  Z channel  data  to  estimate  both  temperature 
and  water  vapor  parameters  has  the  advantage  of  automatically  incorpor- 
ating the  temperature  level  information  for  the  extraction  of  water 
vapor  data  as  proposed  in  this  report,  plus  automatically  providing  the 
water  vapor  corrections  for  the  temperature  sounder  channels. 
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APPENDIX  A 

SOFTWARE  DESCRIPTION 


Loop  Through  Data  Set  Loop  Through  Data  Set 
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Figure  A-l  Initial  Transmittance  Processing 


Loop  Through  Transmittances  Loop  Through  Channels 


t isuro  A -2  Filter  Modified  Transmittance  Data  Set  Generation  Scheme 
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Figure  A-3  Simulated  Radiance  Calculations 
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D MATRICES 
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D MATRICES  (CONTINUED) 
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